Engineering exotic phases for topologically-protected quantum computation by 

emulating quantum dimer models 
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We use a nonperturbative extended contractor renormalization (ENCORE) method for engineer- 
ing quantum devices for the implementation of topologically protected quantum bits described by 
an effective quantum dimer model on the triangular lattice. By tuning the couplings of the device, 
topological protection might be achieved if the ratio between effective two-dimer interactions and 
flip amplitudes lies in the liquid phase of the phase diagram of the quantum dimer model. For a 
proposal based on a quantum Josephson junction array [L. B. loffe et al, Nature (London) 415, 503 
(2002)] our results show that optimal operational temperatures below 1 mK can only be obtained if 
extra interactions and dimer flips, which are not present in the standard quantum dimer model and 
involve three or four dimers, are included. It is unclear if these extra terms in the quantum dimer 
Hamiltonian destroy the liquid phase needed for quantum computation. Minimizing the effects of 
multi-dimer terms would require energy scales in the nano-Kelvin regime. An alternative imple- 
mentation based on cold atomic or molecular gases loaded into optical lattices is also discussed, and 
it is shown that the small energy scales involved — implying long operational times — make such a 
device impractical. Given the many orders of magnitude between bare couplings in devices, and the 
topological gap, the realization of topological phases in quantum devices requires careful engineering 
and large bare interaction scales. 

PACS numbers: 03.67.Pp,74.81.Fa,75.10. Jm 
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I. INTRODUCTION 

Systems characterized by topological quantum order 
(TQO) have a degenerate ground state, which is not 
associated with any broken symmetry, i.e., the differ- 
ent degenerate ground states are indistinguishable un- 
der the action of any local operator^ Instead, they can 
only be distinguished via global operators intimately re- 
lated to their topological properties. TQO does not fit 
into Landau's paradigm for ordered phases of matter^ 
which makes it intrinsically interesting. Furthermore, 
this robustness against local perturbations characteris- 
tic of systems exhibiting TQO can be used to implement 
a fault-tolerant quantum computer^ 

Within this approach, robust storage devices for 
quantum states ("protected memory qubits") can 
be built from Abclian topological quantum states, 
whereas topologically-protected computations ("pro- 
tected gates") can be implemented using non-Abelian 
states^ Given the enormous challenges involved in build- 
ing conventional quantum computers caused by the de- 
coherence inherent to quantum-mechanical systems, the 
alternative approach exploiting topological order has at- 
tracted considerable interest recently because local oper- 
ators (i.e., noise) do not disturb the topological phase. 

One promising class of systems exhibiting TQO are 
fractional quantum Hall systems with filling factors v = 
5/2 and v = 12/5 which are conjectured to exhibit 
non-Abelian anyonic excitations^ Unfortunately, despite 
some evidence,— the existence of anyons in these sys- 
tems remains to be confirmed experimentally. On the 
other hand, a number of interesting lattice models is 



known to exhibit TQO. Among these are quantum dimer 
models (QDM) (Refs. HQIlIi) spin models, and Hub- 
bard models with generalized interactions defined on 
Kagome lattices ) 10 ' 11 ^ 2 ' 13 toricf 3 - and color— codes, as 
well as Kitaev's honeycomb anisotropic spin modelJ^ In 
general, these lattice models incorporate unrealistic ele- 
ments such as artificially-constrained degrees of freedom 
or nontrivial interactions and thus experimental realiza- 
tions remain elusive. Therefore, we are interested in en- 
gineering topologically-ordered phases by emulating lat- 
tice models using highly manipulablc quantum tool-boxes, 
such as Josephson junction arrays^ and cold atomic^ or 
molecular— ^ gases loaded into optical lattices. How- 
ever, as promising as these approaches might seem, the 
challenges imposed to the engineering of such emula- 
tors are huge, requiring special attention to the design 
of such devices and a careful analysis of the involved en- 
ergy scales as well as the possible existence of extra terms 
in the emulated Hamiltonian. 

Having these issues in mind, we use a nonperturba- 
tive algorithm, extended contractor renormalization (EN- 
CORE) (Ref. [190 an extended version of the Contractor 
Renormalization (CORE) tcchniqu o 20 ' 21 to design exotic 
phases to build topological quantum computers as well as 
to propose controllable experiments to investigate TQO. 
We consider an emulator for the QDM on the trian- 
gular lattice based on an array of quantum Josephson 
junctions^ This system is a good candidate for the imple- 
mentation of a topologically-protected qubit for two rea- 
sons: First, quantum dimer models are among the best 
understood systems exhibiting TQO and the presence of 
a topological phase has been unequivocally established in 
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a number of studies^ ' 22 ' 23 ' 24 Second, the manipulation 
of Josephson junctions is an experimentally mature field 
where an exquisite degree of control has been achieved. 
We are able to derive the couplings in the effective model 
describing the low-energy physics in the array in an unbi- 
ased way (the only limitations being caused by the finite 
sizes of the clusters analyzed). Our final conclusion is 
that although the approach of Ref. d seems promising 
based on simple estimates, the energy scales obtained in 
the full analysis are too low to make this approach feasi- 
ble. 

In addition, we also discuss, by means of a per- 
turbative analysis, an implementation based on cold 
atomic/molecular gases loaded into a Kagome-shaped op- 
tical lattice and encounter similar problems of too low 
energy scales and too long time scales. 

II. DEVICES FOR EMULATING QUANTUM 
DIMER MODELS 

A. Quantum Dimer Model on a Triangular Lattice 

The QDM has first been introduced by Rokhsar 
and Kivelson£ in the context of the resonating valence 
bond (RVB) scenario for cuprate superconductors^ The 
square lattice version of this model only displays valence 
bond crystal phases, with the notable exception of a sin- 
gle point at which the correlations decay algebraically 
with distance and the ground state splits into many topo- 
logical sectors^ Its triangular lattice version, first an- 
alyzed by Moessner and Sondhi,— has a gapped liquid 
phase with exponentially decaying correlations extend- 
ing through a finite range of the model parameters. 

The triangular-lattice QDM is given by the following 
Hamiltonian 

Tl = He? + + , (1) 

with 

W Z7 = -^[|-)(//| + |//)<-|] 

(2) 

+^[|-)(_-| + |//}(//|] 
a 

and similar definitions for Tt^ and H.q . Parallel dimers 
sitting on the same rhombus (henceforth we refer to such 
configurations as flippable rhombi) flip with amplitude t 
and interact with each other via a potential strength v; 
the sum runs over all the rhombi with a given orientation. 

Despite its apparent simplicity, the phase diagram of 
the QDM on the triangular lattice is rich, comprising 
different crystalline phases . 7 ' 8 ' 22 ! 23 ' 24 Here we are only 
interested in the quantum liquid phase, which is stabi- 
lized in the range 0.82 < v/t < 1^ with exponentially- 
decaying correlations between dimers and a gap A ~ O.li 
against excitations^ In this phase the system's ground 
state is degenerate: twofold degeneracy on a cylindrical 



geometry, fourfold on a torus (full periodic boundary con- 
ditions). The topological sector to which a given dimer 
configuration belongs can be determined via the parity of 
the dimer count along an arbitrarily chosen reference line 
(see Fig. [l| , a property which can be used to build a two- 
level system for a topologically protected quantum bit 4 
Note that the topologically-ordcred phase of the QDM 
is also stable towards the presence of disorder^ a partic- 
ularly useful feature since the presence of imperfections 
would be unavoidable in any putative engineered device. 



B. Emulator Based on Josephson Junctions 

The emulation of the quantum dimer model on the tri- 
angular lattice can be achieved by using Josephson junc- 
tion arrays. Ioffe et al& introduced two different Joseph- 
son junction array emulators for the QDM. In this work 
we discuss the implementation defined on the Kagomc 
lattice only, since it has a smaller number of supercon- 
ducting islands attached to each site of the underlying 
triangular lattice and thus is more amenable to numer- 
ical studies. However, our main conclusions are imme- 
diately extended to the alternative implementation on a 
decorated triangular lattice. 

The proposed emulator is built from an array of X- 
shaped superconducting islands structured as a kagomc 
lattice, see Fig. Q] (thick black lines). Each X-shapcd 
island is coupled to its four neighboring islands by a ca- 
pacitance Ch and a Josephson current Jh- Inside every 
hexagon of the kagome lattice a star-shaped island (thin 
black lines in Fig. [1]) is placed which couples only ca- 
pacitively to the X-shaped islands via the capacitance 
Cj. The ground capacitance of an X-shaped island is Cx, 
whereas the ground capacitance of a star-shaped island 
is C*. The energies associated with these couplings are 
given by 

where e denotes the elementary charge, and E J = HJ/2e. 
We set H=2e=l. 

One dimer in this array is equivalent to a Cooper pair 
sitting on one of the six X-shaped islands surrounding 
a given star-shaped island, each one corresponding to 
one of the links of the underlying triangular lattice (see 
Fig. [H shaded thick lines). By applying a global bias to 
the array, only half of a Cooper pair is made available 
per star-shaped island and, in order to impose the dimer 
hard-core condition and emulate the QDM Hamiltonian, 
we need to tune the different capacitances and Josephson 
currents. In order to guarantee that no hexagon can par- 
ticipate in the formation of more than one dimer (repre- 
sented by ellipses in Fig. [1]) , we assign a large value to the 
capacitance C\ such that there is a strong repulsion be- 
tween Cooper pairs placed around the same star-shaped 
island. 
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FIG. 1: (Color online) Array of Josephson junctions used to 
emulate the quantum dimer model on the triangular lattice. 
The array is formed by X-shaped superconducting islands 
(thick black lines), which form a Kagome lattice and normal- 
state star-shaped islands (thin black lines) placed at the cen- 
ter of every hexagon of the Kagome lattice. The shaded lines 
are guides to the eye to emphasize the underlying triangu- 
lar lattice of the effective QDM. Cooper pairs hop between 
nearest-neighbors X-shaped islands with an amplitude given 
by the Josephson current J^. A large ratio between the ca- 
pacitances Ci and Ch defines a sizable on-hexagon repulsion 
Z?hcx to emulate the hard-core dimer constraint. The dimers 
are represented by ellipses sitting on one of the six links of a 
given star-shaped island. The parity of the dimer count along 
a reference line f2 (dotted line) is invariant under the dimer 
flips in the Hamiltonian [Eq. |T}]. 



The energetic cost for placing two bosons around the 
same hexagon -Ehex defines the basic energy scale of the 
array. It should not be confused with the on-site re- 
pulsion between two Cooper pairs sitting on the same 
X-shaped island. The parity of the dimer count along 
the reference line (dotted line in Fig. [1]) is invariant 
under dimer flips (local perturbations) in the Hamilto- 
nian [see Eq. {TJ)] and allows for the determination of the 
topological sectors necessary to define a qubit state. 



Emulator Based on Cold Atomic/Molecular 

Gases 



We also consider an implementation of the QDM based 
on cold atomic/molecular gases loaded into a Kagome 
optical lattice, which can in principle be created by using 
three laser beams, 26 with the; following Hamiltonian: 



ft = -W Yl _ n «') + ~T~Y1 n o{no - 1) 



o 



-J$>l&i+&t& i ). 



(4) 



Here rii = b\bi is the bosonic number operator at the 
site i of the Kagome lattice, U is a repulsion between 
two bosons sitting on the same site and J is the hopping 
amplitude between nearest-neighbors sites in the 

Kagome lattice. Eq is the energy required for placing 
two bosons on different sites around the same hexagon 
in the Kagome lattice and enforces the hard-core dimer 
condition, no is the number of bosons sitting around 
a given hexagon. Due to the short-ranged interactions 
between cold atomic gases, the engineering of interaction 
terms as in Eq. §Z§ would likely be a highly nontrivial 
task. One possible solution to this problem is to use polar 
molecule a 17 ' 18 whose permanent dipole moment permits 
long-range interactions. 



III. EFFECTIVE HAMILTONIANS FROM THE 
ENCORE METHOD 

The CORE method was originally introduced by Morn- 
ingstar and Weinstei n 20 ! 21 and since then has been suc- 
cessfully applied to different problems in strongly corre- 
lated systemSi 27 ' 28 i 29 ' 30 i 31 ' 32 i 33 For our application we use 
an extended version, ENCORE, suitable for constrained 
models, such as the quantum dimer models 

The fundamental idea behind CORE and ENCORE 
is to derive an effective model describing the low-energy 
physics of a lattice Hamiltonian by reducing the num- 
ber of degrees of freedom. The usefulness of the method 
relies on a fast decay of the effective interactions for the 
specific effective model, something which needs to be ver- 
ified for each case. A large amount of physical intuition 
is required to obtain physically sound results, which is 
one reason why CORE has not found a more widespread 
use to date. 

The effective Hamiltonian obtained with ENCORE 
generally includes arbitrarily-ranged terms. Large cou- 
plings associated with long-range terms indicate that 
the restricted subspace does not accurately describe the 
low-energy behavior of the original model. However, if 
we are interested in engineering an emulation of a cer- 
tain Hamiltonian, the aforementioned problems are irrel- 
evant because in this case the effective model and the 
restricted Hilbert space are known a priori. If the EN- 
CORE method fails we simply conclude that emulation 
is not possible. 

The breakdown of the mapping is also signaled by the 
appearance of "intruder" states in the low-lying spec- 
trum. These are states with negligible overlap with any 
of the desired low-energy states. Since both aforemen- 
tioned effects are correlated^ we avoid the adoption of 
an arbitrarily-defined threshold value for the long-range 
interactions and we define the breakdown of the mapping 
onto a QDM as the point where a first intruder state ap- 
pears in the device's low-energy spectrum. 

Since our primary goal in the present paper is to verify 
the feasibility of a fault-tolerant quantum bit engineered 
from a system with a topologically-ordercd phase, the 
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device's parameters must be tuned in order to ensure that 
the emulated model has couplings known to correspond 
to a quantum dimer liquid phase. In addition, a careful 
analysis of the involved energy scales is necessary in order 
to avoid technological limitations. 



IV. EMULATING QUANTUM DIMER MODELS 
USING JOSEPHSON JUNCTION ARRAYS 



The array of Josephson junctions discussed in Sec. Ill Bl 
can be described by the following generalized Bosc- 
Hubbard Hamiltonian 



H 



-E 

2 ^ 



n j c i,k nk ~ Jh 



E 



(5) 



The positions of the X-shaped islands in the array are 
denoted by the indices j and k. (J, k) represent nearest 
neighbor (NN) sites in the Kagome lattice, rij = b^bj 
is the bosonic occupation number at site fj, Jh is the 
Josephson current between two X-shaped islands. 

C _1 is obtained by numerically inverting the capaci- 
tance matrix C of the array. The matrix elements con- 
necting two X-shaped islands in this matrix are given by 



Cj.k — [Cx + MjC; + VjCh\5 r -^ T r- k + GhSf j>r - k 



(G) 



where f connects NN sites in the Kagome lattice, fx is the 
number of hexagons a given X-shaped island joins [/j = 2 
for full periodic boundary conditions (PBC)] and v is its 
number of NN [v = 4 for PBC]. 

The normal-state star-shaped islands are only capaci- 
tively connected to the X-shaped islands and their only 
role is to set up the interactions in the Hamiltonian. The 
inverse C _1 appearing in the Hamiltonian is sensitive 
also to these interactions, specified in the following. Star- 
shaped islands sitting on the sites R a and Rp of the un- 
derlying triangular lattice contribute with, 



(I) 



and the elements connecting X- and star-shaped islands 
are 



Cj,a — Ci^rj ,ra+s*: 



(8) 



where s are the vectors connecting a star-shaped island 
to the X-shaped islands surrounding it. 

The energy -Ehcx to place two dimers on a hexagon 
can be obtained from certain matrix elements of C _1 . 
Quantum fluctuations due to the Josephson coupling Jh 
reduce this bare value and we thus include them in second 
order in perturbation theory in our discussions below. 
To ensure that we are allowed to restrict the calculations 
to hard-core bosons, we have verified that the on-site 
repulsion is larger than Jh by a factor of at least 50 for 
all sets of couplings in the array. 
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FIG. 2: (Color online) Open-boundary clusters studied: (a) 
N x2 (N = 3 in the figure) hexagon ladders; (b) ten-hexagon 
cluster; (c) - (e) special clusters with four, six and eight 
hexagons, accommodating the lowest-order flip (represented 
by the associated transition graphs where full dimers flip to 
open ones) involving two, three, and four dimers, respectively. 



Our results arc obtained by analyzing the small open- 
boundary clusters depicted in Fig. [5J ladder-like clus- 
ters with N x 2 hexagons (N = 3, 4, and 5), a ten- 
hexagon cluster (from which most results have been ob- 
tained), and three special clusters with four, six, and 
eight hexagons which accommodate only two distinct 
dimer configurations each. 

Finally we need to take into account experimental 
limitations. The smallest values for capacitances be- 
tween two superconducting islands obtainable with cur- 
rent technologies are such that E c = 1/2C ~ IK [see 
Eq. ([3])], higher values of E c can be obtained for ground 
capacitances. We set the smallest junction capacitance 
Ch = 0.5, such that E£ = 1 to set the energy scale, and 
we restrict our analysis to values of C\ > Ch throughout 
the rest of this paper. In this way, assuming a value of 
E^ rj IK, a priori taking into account current techno- 
logical limits, all energies are fortuitously directly given 
in Kelvin. 



A. Dimer Flips 

The simplest dimer flip involves two parallel dimers on 
the same rhombus of the triangular lattice, as illustrated 
in Fig. [He). It involves the creation of a virtual state 
in which one hexagon is doubly occupied, occurring with 
an amplitude given in second-order perturbation theory 
by t « Jh/^hcx- Using ENCORE, we now analyze the 
amplitude associated with this dimer move for the set of 
capacitances studied in Ref. d: 



C» = 10, C x 
Q = 2.0, Ch = 



= 10, 
0.5. 



(9) 
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FIG. 3: (Color online) (a) Two-dimer flip amplitude t cal- 
culated with ENCORE versus the Josephson current Jh for 
the clusters depicted in Fig. [2] The parameters of Ref. @ are 
used: C» = C x = 10, G = 2, and C h = 0.5. Solid lines 
are only guides to the eye. The vertical dashed line indicates 
the point where the mapping onto a QDM breaks down (cf. 
Sec. IIIip . Second-order perturbative results obtained by nu- 
merically calculating I?hex are indicated by the thick black 
curve. The dotted curve corresponds to the results obtained 
by using the expression for Shex derived in Ref. (b) Ampli- 
tudes for each of the five nonequivalent two-dimer flips in the 
ten- hexagon cluster (dashed curves, see main text), compared 
to their average (downward triangles). 



In Fig. [Ha) we show the results for the clusters de- 
picted in Figs. [Ha)-{5Jc). The absence of substantial 
finite-size effects confirms that the two-dimer flip is a 
local process. The different clusters differ only after 
the point where the mapping onto a QDM breaks down 
due to the appearance of intruder states (see Sec. IIII[) . 
The results from the ten-hexagon cluster are obtained 
from an average between the amplitudes for five possible 
nonequivalent two-dimer flips occurring within slightly 
different "dimer environments," the different ways the 
dimers not participating in the flip arc arranged in the 
cluster J£ Again, the amplitudes for these individual pro- 




C h (K ) 

FIG. 4: Numerical second-order perturbative results for the 
on-hexagon repulsion -Ehcx for the ten-hexagon cluster as a 
function of the capacitance Ch for C, = Cx = 10 and G = 
2. The horizontal dotted line is _B hex ~ 0.2(Ci/C x ) 2 K? , as 
obtained in Ref. M. 



cesses only deviate slightly from their average until the 
point where the dimer picture breaks down, as shown in 

Fig.ETb). 

In Fig. [3f a) we also compare to second-order pertur- 
bative data for t, obtained from t = J^/£"hex by numer- 
ically calculating i?hex for the ten-hexagon cluster (solid 
line; see Fig. [5]), as well as by using the approximation 
£ hGX - 0.2(d/C x ) 2 EC of Ref. H The discrepancy be- 
tween the two estimates clearly illustrates the nontrivial 
dependence of -Ehcx on the set of capacitances adopted for 
the array. Note that -Ehcx vanishes when Ch = C\. Us- 
ing the accurate estimate for i?hcx (solid line in Fig. Ql, 
we find reasonable agreement with the ENCORE results, 
which motivates us to use the perturbative results to 
guide our optimizations below. 

Additional flips involving three and four dimers occur- 
ring within third and fourth order in the small parameter 
Jh with amplitudes t% and t± respectively, are depicted 
in Figs.[2td) and[2^e). Together with the two-dimer flip, 
these processes are special because they are the lowest- 
order possible dimer moves in the array: comprising m 
dimers, they appear as mth order processes in Jh- All 
other flip terms are strongly suppressed in the limit of 
small currents Jh- Throughout this paper, we denote the 
sum of the absolute values of the amplitudes for all other 
dimer moves by S, and we use this quantity to gauge the 
validity of the mapping onto a QDM. 

Of particular interest is the flip involving three dimers 
in a triangular configuration depicted in Fig.OJd), which 
has also been found in a recent mean-field mapping by 
Vernay et al^ of a spin-orbital model for the compound 
LiNi02 onto a QDM. They found that this extra dimer 
move considerably to extend the liquid phase, allowing 
for "extra room" in trying to optimize the couplings in 
the array. One is tempted to conclude that the four-dimer 
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FIG. 5: (Color online) (a) Total Coulomb energy versus num- 
ber of flippable rhombi (parallel dimers which can flip) for a 
4x4 PBC cluster in the limit J h = and with C* = Cx = 10, 
d = 2. and Ch = 0.5. Deviations from the dominant linear 
behavior signal the presence of extra interaction terms beyond 
the one between parallel dimers sitting on the same rhombus 
with strength v. (b) Interactions present in the limit Jh = 0. 
In addition to the dominant rhombus term v, two further 
interaction terms involving three dimers with amplitudes u\ 
and U2 have to be included. The deviations from the main 
contribution in panel (a) can be explained by the presence of 
these extra terms. 



ber of flippable rhombi, confirming that the interaction 
between parallel dimers sitting on the same rhombus v 
is the dominant diagonal term in the emulated dimcr 
Hamiltonian. The deviations from the linear behavior in 
Fig. [UJa) indicate that other interaction terms are also 
present. These extra contributions cannot be explained 
by pairwise interactions, but all such deviations are en- 
tirely described if we take into account three-dimcr inter- 
actions with strengths u\ and u-i [as shown in Fig. EJb)]. 

Understanding the dependence of the interaction 
strengths on the capacitances of the array is essential 
for the optimization of the energy scales in the emulated 
Hamiltonian, and for the minimization of the couplings 
associated with the three-dimer interactions, the effect of 
which arc unknown. Figure [5] shows the dependence of 
the couplings v, ui, and u-i on the capacitances^ on the 
same 4x4 cluster with PBC for Jh = 0. The interac- 
tion term v( Jh = 0) peaks for small values of C\ and Cx, 
and iti (Jh = 0) and U2(Jh = 0) are appreciable over a 
significant regime of the parameter space. Decreasing C* 
leads to larger values of «i(Jh = 0) and 1/2 (Jh = 0) (not 
shown) but has no significant effects on u(Jh = 0). 

To include the effects of nonzero Josephson current 
Jh on the interactions, we study the ten-hexagon clus- 
ter shown in Fig. [^b) using ENCORE. Due to the small 
size of the cluster, only four-dimer configurations, out of 
a total of 14 in this cluster, give nonequivalent diagonal 
contributions in the effective dimer Hamiltonian. Bound- 
ary interactions due to open boundary conditions further 
complicate the estimates. Without going to larger clus- 
ters (which is not possible with current computational 
resources) , we cannot determine the individual terms but 
only the following combined quantity: 

v' =v- («! - u 2 ) • (10) 

which is nevertheless a good estimate for v as we can see 
by analyzing the 4x4 cluster. As shown in Fig.^d), v' 
underestimates v by typically only 20% for Jh = in the 
region where we have a valid mapping onto the QDM, and 
these corrections do not change the conclusions drawn 
below 1^ 



term £4 has similar effects, and it will be of interest to 
confirm this numerically. 



B. Dimer Interactions 

Dimcr-dimer interactions have a nontrivial dependence 
on the particular choice of capacitances in the array and 
must be tuned in order for the ratio v/t to lie in the 
liquid phase. To investigate them, we first analyze a 
4x4 cluster with PBC in the limit of zero Josephson 
current ( Jh = 0) by inverting the capacitance matrix on 
this cluster. In Fig. [5£a) we show the Coulomb energy 
as a function of the number of flippable rhombi in each 
configuration. The total energy scales well with the num- 



C. Tuning the Ratio v/t 

Having discussed the dimer flips and interactions in the 
array of Josephson junctions, we now analyze the feasi- 
bility of the emulation of a topological phase by adjusting 
the couplings and required energy scales in the emulated 
QDM to achieve an effective model with couplings in the 
desired range 0.82 <v/t<l. 

In Fig. [7| we show the ratio v'/t as a function of the 
Josephson current Jh for the set C* = Cx = 10, C\ = 2 
and Ch = 0.5 as proposed in Ref. d (|36h . As seen in 
Fig. the topological phase corresponding to v'/t < 1 
(marked by dashed horizontal lines) is not reached before 
the breakdown of the mapping to the QDM. This break- 
down, marked by a vertical dashed line is seen both in 
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FIG. 6: (Color online) Dependence of the amplitudes v, iti and U2 on the capacitances Cx and d, for C» = 1 and Ch = 0.5. 
The results are obtained from the analysis of the 4x4 PBC cluster with Jh = 0. Panels (b) and (c) suggest that minimization 
of ui and U2 is hard to be achieved and panel (d) shows that our results for v' = v — (u\ — U2) (our estimate for v for finite 
values of Jh, see main text) underestimate v by a factor of 40% at most, and typically about 20%. 



the appearance of nondimcr-like states in the low-lying 
spectrum as well as by a drastic increase in the summed 
amplitude for multi-dimer flips (£). We conclude that 
for this value of capacitances no topological phase exists. 

In order to stabilize the mapping onto a QDM for 
larger values of the Josephson current, we explore al- 
ternative sets of capacitances leading to larger values for 
£ticx, more strongly suppressing non-hardcore-dimer con- 
figurations. The goal is to avoid the breakdown of the 
mapping before the target values for v/t — where the sys- 
tem is in the topological liquid phase — is reached. The 
dependence of -Ehex on the capacitances C\ and Cx for 
two arbitrarily chosen values of C* , obtained from a nu- 
merical second-order perturbative analysis of the ten- 
hexagon cluster is shown in Figs.E^a) and[Uc). There is 
a region close to the C axis where i?hex peaks, the peak 
value increasing for smaller values of C* . 



Besides maximizing £hcx, we also want to optimize the 
amplitude t since the topological phase appears only for 
temperatures below the gap A ~ O.lt. We use second- 
order numerical perturbative results for the ten-hexagon 
cluster as a guide. For each set of capacitances we calcu- 
late the value of the Josephson current Jh giving the de- 
sired ratio v/t = 1 , and we plot the value of t in Figs. Efb) 
and[8jd). Since within perturbation theory, we are not 
able to determine whether a particular set of parameters 
leads to a valid mapping onto a QDM, we introduce an 
arbitrary cutoff, Jh = Ehcx/2, which, based on our EN- 
CORE results is a generous upper bound. Capacitances 
for which the target ratio 0.82 < v/t < 1 is not reached 
below this bound are discarded and indicated as blank 
regions in the figures. 

These perturbative results show that optimal values 
for t arc obtained in a region corresponding to small val- 
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FIG. 7: (Color online) Our estimate v' = v — (wi— ua) (circles) 
(Ref . l3rjl) for the two-dimer interaction v vs the Josephson cur- 
rent Jh, for C, = Cx = 10, d = 2 and Ch = 0.5, computed 
with ENCORE on the ten-hexagon cluster. Triangles denote 
the sum of the absolute values of the amplitudes (E) associ- 
ated with all possible dimer moves in the ten-hexagon cluster 
beyond the lowest order ones depicted in Fig. [2][c)-[2fe). The 
breakdown of the mapping onto a dimer model, indicated by 
the vertical dashed line, correlates with an abrupt increase 
in E and occurs before the topological phase of the simplest 
QDM for 0.82 <v/t< 1.0 (indicated by the dashed horizontal 
lines) is reached by v'/t. 



for the considered value of C*. Larger values for the 
flip amplitude can be obtained if we choose smaller C*. 
But as shown in Fig. OJb) , only slightly larger values of 
imax ~ 3mK are obtained for values of C* two orders of 
magnitude smaller, suggesting that saturation is rapidly 
reached. Since experimentally we cannot arbitrarily de- 
crease the ground capacitances, we conclude that we can 
estimate an upper-bound for the amplitude t consistent 
with 0.82 < v/t < 1 of only a few milli-Kelvin. These 
optimal values are associated with small Josephson cur- 
rents only slightly larger, Jh < lOmK, typical experimen- 
tal values being close to IK. Returning to the fact that 
v' underestimates v (Sec. IIVB[) we see that this does not 
influence this upper bound estimate. 

Thus, even if we assume that the longer-range terms 
present in the emulated QDM do not destroy the topo- 
logical liquid phase and that we can still estimate the 
topological gap as A ~ O.li, we can expect that the op- 
erational temperatures for the putative quantum bit is in 
the micro-Kelvin regime, clearly far below the limits of 
current technologies. 



D. Extra Flips and Interactions 



ues of Cx and close to the point where the breakdown 
of the mapping onto a QDM occurs before the target ra- 
tio is reached. This behavior can be understood if we 
analyze the dependence of v(Jh = 0) and i?hex on Cx 
and d, shown in Figs. Ufa), Fig. Ufa), and^c). Since 
the ratio v/t decreases monotonically with Jh, larger val- 
ues of u(Jh = 0) have the desirable effect that values 
0.82 < v/t < 1 are reached for larger values of Jh, associ- 
ated with more favorable values for t and A. Decreasing 
C\ at small values of Cx, we can see from Fig. Ufa) that 
progressively larger values of v(Jh = 0) can be obtained. 
However, eventually, the region where i?hex peaks is sur- 
passed and the mapping onto a QDM breaks down before 
the target ratio for v/t is reached. Decreasing C* does 
virtually not affect v( Jh = 0) but considerably increases 
Ehcx, therefore allowing us to obtain even more favorable 
values for t, as shown in Fig. [HJd). 

We have confirmed the validity of this qualitative anal- 
ysis using ENCORE. In Fig. [S] we show v'/t as a function 
of Jh for two different values of C* . For each value of C , 
Cx is chosen such that -Ehcx is close to its maximum and 
optimal values for t can be obtained (see Fig. [8]). We 
terminate each curve when the mapping breaks down, as 
seen by nondimer intruder states in the low-energy spec- 
trum. In Fig. [Ufa) we see that for C» = 1 and Cx = 0.25 
more favorable results for t can be obtained by decreas- 
ing C , but below C ~ 1.5 the mapping onto a QDM 
breaks down before the target ratio 0.82 < v/t < 1 is 
reached. We thus get an upper bound of i max ~ 1.5mK 



So far we have ignored the presence of interactions 
and flips comprising three or more dimers in the effec- 
tive Hamiltonian, although they are most likely relevant 
as suggested by the analysis of the Jh = limit >21 In 
Fig. [TOF a) we show the dependence of the couplings on 
J h for C* = 1, C x = 0.25, d = 2 and C h = 0.5. In par- 
ticular, the amplitudes associated with the unfrustrated 
flips involving three and four dimers [Figs.^d) and^c)] 
are larger than the topological gap A ~ O.lt of the stan- 
dard QDM. Vernay et al££ showed that the three-dimcr 
flip extends the liquid phase. Thus, in order to be able 
to precisely estimate the operational temperatures for the 
emulated qubit, it is necessary to study the effects of the 
inclusion of the extra terms in the QDM. However, even 
in the absence of such a detailed analysis we can conclude 
that the involved technological challenges in reaching the 
sub- milli-Kelvin temperatures required for this device are 
substantial. 

Aiming for a simpler QDM with only two-and three- 
dimer flips [Figs. OJc) and[^d)] requires the suppression 
of higher-order flips. We find that this requires very small 
values of Jh and leads to even smaller values of t. To illus- 
trate this problem, we show in Fig. fTUTb) the dependence 
of all couplings on Jh in the emulated QDM for C* = 100, 
C x = 500, C = 10, and C h = 0.5. With J h = the 
interactions are given by v — 0.879nK, u\/v = 0.006, 
u-i/v = 0.005, and v' = 0.878nK; nano-Kelvin tempera- 
tures are unrealistic in a solid-state device. 
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FIG. 8: (Color online) The on-hexagon repulsion i?hex (left panel) and the dimer flip amplitude t for v/t = 1, obtained by a 
second-order perturbative analysis of the ten-hexagon cluster, as a function of the capacitances Cx and d for C, = 1 and 0.01. 
White regions in panels (a) and (c) indicate capacitances for which -Ehex is not a well defined quantity. In panels (b) and (d) 
white regions correspond to sets of parameters for which the target ratio for v/t is reached for Josephson currents Jh > -Ehex/2, 
an upper bound for the largest current leading to a valid mapping onto a dimer model. 



V. EMULATING QUANTUM DIMER MODELS 
USING COLD ATOMIC/MOLECULAR GASES 

We now turn our attention to the alternative imple- 
mentation based on cold atomic/molecular gases loaded 
into an optical lattice which was presented in Sec. Ill CI 
Since no concrete microscopic proposal is available, we 
restrict ourselves to order of magnitude estimates. 

Flips comprising two dimcrs in this system similarly 
involve the creation of a virtual state with energy Eq and 
thus occur, within second order in J, with amplitude t = 
J 2 /Eq (we can also expect that flips involving three and 
four dimers may play an important role here). Since the 
mapping onto a QDM necessarily breaks down when the 
kinetic energy dominates over the on-hexagon repulsion, 
it follows that an upper bound for the hopping amplitude 



consistent with a dimer picture is given by J ma x ~ Eq/A, 
and the largest obtainable value for the flip amplitude in 
such emulator is thus i max = J 2 lax /Eo w J max /4. 

Preparation of the quantum bit state requires a con- 
trolled mixing of dimer states belonging to different topo- 
logical sectors. This can be achieved by virtually break- 
ing one dimer and creating a virtual particle-hole exci- 
tation, the particle corresponding to a doubly-occupied 
hexagon and the hole to an empty one 4 For a qubit 
with linear dimension corresponding to M hexagons, an 
upper-bound for the mixing amplitude h x is given by 

(j \ M - 1 n\ M - 1 

The largest attainable experimental values for the hop- 
ping amplitude in cold atomic gases loaded into optical 
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FIG. 9: (Color online) Dependence of the estimate for two- 
dimer repulsion v jt on the Josephson current Jh and capac- 
itance Ci for the ten-hexagon cluster obtained with the EN- 
CORE algorithm for (a) C, = 1, Cx = 0.25, and C h = 0.5, 
and (b) C* = 0.01, Cx = 0.05, and C h = 0.5 (Ref.[37|). The 
last data point on the right-hand side corresponds to the value 
of Jh for which the mapping onto a QDM breaks down. The 
dashed lines mark the spin-liquid phase with TQO for the 
standard QDM. 



lattices are close to 1 kHz, smaller values being expected 
for more massive molecules. Thus, even on a rather small 
lattice comprised of the 10 x 10 hexagons (M Kt 10), 
we can conclude that the time-scale involved in a sin- 
gle qubit manipulation is of the order of minutes, much 
longer than typical coherence times in cold atomic gases 
in optical lattices. 



VI. CONCLUSIONS 

We have studied proposals to emulate a triangular lat- 
tice quantum dimer model (QDM). A realistic emula- 
tion of the QDM would allow the implementation of a 
fault-tolerant quantum bit, allowing us to circumvent the 
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FIG. 10: (Color online) Flip amplitudes and interactions ob- 
tained by ENCORE for (a) C* = 1, Cx = 0.25, C = 2, and 
C h = 0.5 and (b) C* = 100, Cx = 500, C = 10 and C h = 0.5. 
In the upper panel only, the rightmost data point indicates 
the breakdown of the mapping onto a QDM . Dashed hor- 
izontal lines indicate the range of parameters corresponding 
to the topological phase of the standard QDM. Amplitudes 
for the special flips involving three (t^) and four (£4) dimers 
are also shown. For {3, results are obtained from the analysis 
of the ten (filled triangles) and six (empty triangles) hexagon 
clusters; £4 is calculated from the eight-hexagon cluster. In 
the limit Jh = the strengths of dimer interactions are: (a) 
v = 525/iK (v' = 390^K), ui/v = 0.512, and u 2 /v = 0.253 
and (b) v = 0.879nK (v' = 0.878nK), ui/v = 0.006, and 
U2/1) — 0.005. Values for t corresponding to the target ratio 
0.82 < v/t < 1.0 are t ~ lmK for the set in panel (a) and 
t ~ 2nK for the set in panel (b). 



problem of decoherencc which plagues more conventional 
proposals for achieving quantum computation. 

The Josephson junction emulator by Ioffe et was 
studied numerically using the ENCORE method. Our 
results showed that the largest attainable values for the 
two-dimcr flip amplitude t are a few milli-Kelvin and re- 
quire very small Josephson currents. Since a device based 
on such an array would only be operational at temper- 
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atures considerably below the topological gap A ~ O.lt, 
implementation of a topologically-protected quantum bit 
with the considered array is beyond the present day tech- 
nology. 

The alternative array introduced by Ioffe and 
collaborators^ comprised of Y-shaped superconducting 
islands forming a decorated triangular lattice would lead 
to even lower values for the flip amplitudes since dimers in 
this implementation correspond to a resonating Cooper 
pair and dimer flips involve the tunneling of this pair 
through a weaker link with much smaller Josephson cur- 
rent. Similar challenges with too low energy scales and 
too long time scales are also faced by implementations 
using cold atomic gases. 

Our results illustrate the challenges involved in the de- 
sign of emulators for exotic phases. The fundamental 
reason for these difficulties resides in the fact that topo- 
logical quantum order is a low-temperature feature, since 
the system's local degrees of freedom must be highly en- 
tangled over long distances, of order of the system's size, 
for topological order to emerge. Since emulation of the 



relevant models is obtained in the low-energy limit of 
the proposed quantum device, we face the challenge that 
extremely low temperatures are required. Thus, the ap- 
proach of emulating topologically ordered states for per- 
forming fault-tolerant quantum computation might only 
be a successful one if we can devise emulators based on 
much stronger bare electronic interactions, and a detailed 
analysis of engineering limits is required. 
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